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La presents invention concerne un procede at un dispositif 
d'evaluation de la qualite des images codees ainsi que Tutilisation de tels 
procede et dispositif. 

Dans les systemes de codage de sequences video numeriques, il 
5 est connu d'estimer la qualite des images en sortie de processus par 
comparaison avec Timage originate en utilisant le rapport signal a bruit. 

Ce rapport est en general denomme PSNR (Peak Signal to Noise 
Ratio) et obtenu en sommant les differences quadratlques des pixels de 
rimage finale et de Timage originate. 
10 Toutefois. cette mesure ne prend pas en compte les caracteristiques 

psychovisuelles de la vision humaine (HVS : Human Visual System), En 
effet, Toeil humain est plus sensible a certaines frequences spatiales et sa 
perception des details des objets est fortement liee a leur mouvement relatif 
et egalement aux phenomenes de luminosite et de contraste. 
15 Ainsi une sequence, qui peut. selon le processus d'estimation de 

qualite traditionnelle. apparaitre comme resultant d'un bon codage et etre 
supposee avoir une bonne qualite d'images, ne sera pas perdue de cette 
iagon par Tobservateur compte tenu de ses caracteristiques psychovisueltes. 

L'objet de invention est done de permettre une appreciation de la 
20 qualite des images qui se rapproche autant que possible de la perception 
qu'en aura Tobservateur. 

L'influence du facteur humain (HVS) est partietlement prise en 
compte dans les codeurs traditionnels de type MPEG2 dans la perception 
des frequences spatiales des images decodees par utilisation d'une matrice 
25 de ponderation sur les hautes frequences des blocs d'images 8 par 8. mais 
cette prise en compte n'est absolument pas enectuee aans la perception des 
details des objets en mouvement. 

La plupart des travaux mesurant les disparites entre deux images 
sont essentiellement bases sur une analyse statique des d6fauts ne prenant 



pas en compte le mouvement, ou encore sur une analyse des frequences 
spatiales. 

Par ailleurs. les quelques travaux anterieurs prenant en 
consideration I'aspect tenriporel, mais ne prenant pas en compte le 
mouvement proprement dit, utilisent la difference (distorsion) entre un 
macrobloc courant et le macrobloc a la meme position dans Timage 
precedente. 

Aucun de ces travaux ne considere I'influence de la vision humaine 
avec les problemes de mouvement dans Timage. 

Le but de rinvention est de pallier les inconvenients de Tart anterieur. 

Ce but est atteint par le fait que le precede d'evaluation de la qualite 
dimages comporte : 

a) une etape de traitement du signal representatif de Timage pour 
obtenir un signal traite, 

b) une etape de construction a partir du signal representatif de 
rimage codee, d'un signal representatif de I'image champ de mouvement a 
partir de la sequence source, 

c) une etape de constitution d'un signal representatif de la 
segmentation du champ de mouvement et de stockage des pixels d'images 
representatifs de chaque region ayant un champ de mouvement different a 
une adresse definie par rapport aux vecteurs Vitesse estimes dans Tetape de 
construction du champ de mouvement permettant de determiner les pixels 
ayant des vecteurs Vitesse differents, 

d) une etape de determination ou de calcul d'un filtre humain 
psychovisuel a appliquer en fonction de la vitesse estimee de la region, 

e) une etape de filtrage du signal traite, et 

f) une etape de construction de la carte des disparites entre les 
signaux representatifs de Timage obtenus apres I'etape de filtrage et les 
signaux representatifs de I'image decodee obtenus apres Tetape de filtrage. 




Selon une autre particularite, le precede d'evaluation de la qualite 
d'images comporte une etape consistant a appliquer chacune des etapes a 
rimage source et a Timage d6codee. 

Selon une autre particularite, le precede d'6valuation de la quality 
5 d'images codees comporte une etape de decomposition frequentielle des 

images (FFT. sous-bande ) prec6dant Tetape de filtrage constitute d'une 

ponderation par un coefficient d6duit de courbes prenant en compte la 
Vitesse estimee et la bande de frequence consid§ree. pour tenir compte de 
influence relative de la vitesse et de la frequence spatiale sur la perception 

10 des images animees. 

Selon une autre particularite. I'etape de filtrage psychovisuel 
s'applique des matrices representatives des differences inter-pyramides 
entre les pyramides de Laplace des images sources traitees et celles des 
images decodees traitees apres ponderation par, d'une part Tinfluence locale 

15 representative de la frequence du pixel concerne et d'autre part, un 
coefficient de filtrage deduit de courbes de filtrage prenant en compte la 
Vitesse estimee et la bande de frequence correspondent au niveau de la 
pyramide de Laplace a laquelle le pixel appartient dans une pyramide 
muitiresolution obtenue par constmction d'une pyramide a partir de I'image 

20 de chaque r6gion de vitesse differente. 

Selon une autre particularite, les courbes de filtrage psychovisuel 
sont soit constituees par une succession de courbes arrangees sous forme 
de base de donnees et memoris§es dans le systeme, et 6ventuellement 
interpolation a partir de ces courbes, soit obtenues par une representation 

25 analytique mise en oeuvre par des moyens de calcul permettant de calculer 
chaque courbe. 

Selon une autre particularite, Tetape de construction de la carte des 
disparites s'effectue par recomposition des pyramides multiresoiutions 
filtrees obtenues a retape precedente. 



Selon une autre particularity, r§tape de traitement de I'image 
comporte une etape de decomposition des images sources et d§codees en 
pyramide de Laplace de n niveaux et une etape de construction de la 
difference inter-pyramide. 
5 Selon une autre particularity, la Vitesse ou valeur locale du 

mouvement est obtenue par une construction eventuelle de filtres puis 
application du filtre construit ou par application d'un filtre median. 

Selon une autre particularite, le precede comporte une etape de 
pr§correction des images en effectuant une correction Gamma et une 
10 correction par la loi de Weber. ^ 
Selon une autre particularity, la correction Gamma s'effectue par un 
dispositif mettant en ceuvre la fomnule suivante : 

_y^_Ks-V^ " 

avecV=kaE^ 

15 dans laquelle y est la luminance, V la tension de luminance, E I'ydairement 
de I'image analysee y est une exposant voisin de 2,2 pour les tubes images 
en noir et blanc et ya une valeur de 0,45 admise couramment pour la 
television couleur. 

Selon une autre particularity, le filtrage est obtenu par construction 
20 du filtre psychovisuel correspondant a la vitesse estimye d partir d'une base ^ 
de donnyes de filtres et interpolation entre les deux filtres correspondant aux 
rygions les plus proches de la rygion dont la vitesse a yty estimye. 

Selon une autre particularity, influence (l„) locale relative du pixel Pi 
concemy est obtenue par calcul d'une valeur En reprysentant la puissance 
25 q'^"** de la diffyrence niveau d niveau inter-pyramides entre les pyramides 
sources et pyramides decodees de meme niveau du pixel concemy. 

Selon une autre particularity, le calcul de In est effectue par 
utilisation de la formule suivante 



avec E„ = (DiffnCPij))" . 
m(EK)= Ek si Ek > S 
et m(Ek) = S si Ek < S 

avec par exemple S = 0,5% (valeur maximale possible de Ek) 

Selon une autre particularlte, le filtrage comporte un filtrage 

directionnel des images selon une direction d6termin6e plutSt qu'une autre. 

Selon une autre particularity, la correction Gamma s'efFectue par un 

dispositif de calcul mettant en oeuvre I'^quation suivante : 



10 ^display "^max 



e etant la valeur de niveau de grille du pixel, e^ax ©tan* 'a valeur maximum 
exemple. .25.6_ si le coda g e s'eff e ctue sur 8 bits. Lm ^r 6t ant I'i ntensit6 
correspondant k gfnax ®ri cd/m2. 

Selon une autre particularlte, la loi de Weber est mise en oeuvre par 
15 un dispositif de calcul qui realise la fonction suivante : 



Vo^=%-IoJl-HlOO-f^ 



Selon une autre particularite. le calcul du filtre psychovisuel s'obtient 
20 par un dispositif mettant en oeuvre la formule suivante : 

G(a.v) = [6.1 + 7.3(log(v/3)^] x va^ exp[-2a (v+2)/45.9] 
avec^^=27tfrf=freqtjence^attai^^ 
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Un autre but de Tinvention est de proposer une utilisation du procede 
selon invention. 



Get autre but est atteint par le fait que le precede de I'invention est 
utilise dans un dispositif de codage par retroaction dynamique en fonction 
des disparites psychovisuelles calcul6es par le dispositif de calcul mettant en 
ceuvre le proc6de sur I'un des parametres utilises par le dispositif de codage 
au cours du codage. 

Selon une autre particularity, les disparit§s calcuiees sont 
comparees d un seuil pour modifier les parametres de codage de I'appareil 
de codage jusqu'd ce que le seuil souhaite soit franchi. 

Selon une autre particularity, I'un des parametres de codage est soit 
le pas de quantification, soit la taille des images, soit la forme du groupe 
d'imagesGOP. 

Selon une autre particularite, le precede de I'invention est utilise 
dans I'analyse de rhomog§neite des disparites calcuiees pour agir sur les 
parametres de codage. 

Selon une autre particularity, le precede de I'invention est utilise pour 
modifier les parametres de codage des objets diff§rents d'une image dont le 
codage est oriente objet en fonction d'une disparity souinaitye constante. 

Selon une autre particularity, le precede de I'invention est utilisy pour 
effectuer une ryallocatlon dynamique des dybits alloues. 

Un dernier but de I'invention est de proposer un dispositif mettant en 

ceuvre le procydy. 

Ce but est atteint par le fait que le dispositif d'y valuation comporte : 

- un moyen (1a, 1b) de traitement du signal reprysentatif de I'image 
source (10a) et de I'image dycodye (10b) pour obtenir un signal d'image 
source traitye et d'image decodye traitye, 

- des moyens (2a, 2b) de construction a partir du signal reprysentatif 
de chacune des images, d'un signal representatif de I'estimation du champ 
de mouvement a partir de chacune des images des syquences sources et 
decodees, 




- des moyens (3a. 3b) de constitution d*un signal representatif de la 
segnnentation du champ de mouvement et de stockage des pixels d'images 
representatifs de chaque region Rj ayant un champ de mouvement different 
a une adresse d6finie par rapport aux vecteurs vitesse estimee dans Tetape 
de construction du champ de mouvement permettant de determiner pour 
chacune des images sources et decodees celles ayant des vecteurs vitesse 
differents. 

- un moyen (4,5) de determination ou de calcul d'un filtre humain 
psychovisuel a appliquer en fonction de la vitesse estimee de la region. 

- des moyens (6a, 6b) de filtrage appliques a chacune des images 
sources traitees et decod6es traitees et 

- un moyen (7) de construction de la carte des disparit6s entre les 
signaux representatifs de Timage source traitee obtenus apres Tetape de 
filtrage et les signaux representatifs de Timage decodee traitee obtenus 

15 apres Tetape de filtrage. 

Selon une autre particularite, les moyens de filtrage psychovisuel 
s'appliquent a des matrices representatives des differences inter-pyramides 
calcul6es par des moyens de calcul entre les pyramides de Laplace des 
images sources traitees et celles des images decodees traitees apres 

20 ponderation par, d'une part Tinfluence locale representative de la frequence 
du pixel concerne et. d'autre part un coefficient de filtrage deduit de courbes 
de filtrage memorisees ou calculees et prenant en compte la vitesse estimee 
et la bande de frequence correspondant au niveau de la pyramide de 
Laplace a laquelle le pixel appartient dans une pyramide multiresolution 

25 obtenue par des moyens de construction de cette pyramide multiresolution d 
partir de Timage de chaque region de vitesse diff6rente. 

Selon une autre particularite, les moyens de construction de la carte 
des disparites effectuent une recomposition des pyramides multir§solutions 
filtrees. 





Selon une autre particularite, les moyens de traitement, les moyens 
de constitution, les nnoyens de determination, les moyens de construction, 
les moyens de filtrage sont constitues d'au moins un microprocesseur 
associe a des memoires suffisantes pour contenir les programmes 
permettant de realiser les diff6rents moyens et pour contenir les bases de 
donnees et les informations intermediaires necessaires au calcul et a 
Tobtention de la carte des disparites. 

D'autres particularites et avantages de la presente invention 
apparaitront plus clairement a la lecture de la description ci-apres faite en 
reference aux figures annexees dans lesquelles : 

- la figure 1a repr§sente une vue synoptique des etapes d'une 
premiere variante du precede, 

-la figure 1b represente une representation graphique du resultat des 
etapes de pretraitement, 

- la figure 1c represente une vue simplifiee de la matrice d'image 

source. 

- la figure Id represente une vue simplifiee de la matrice obtenue 
apres segmentation du champ de mouvement, 

- la figure 2 represente une vue synoptique des differentes etapes 
d'une deuxieme variante du precede, 

- la figure 3 represente la famille de courbes de filtrage 
correspondant ^ influence psychovisueile de la vision humaine. courbes qui 
sont memorisees dans une base de donnees pour des vitesses determinees. 

- la figure 4 represente la pyramide multiresolution, 

- la figure 5 represente la pyramide de Laplace. 

Une premiere variante de realisation de Tinvention va etre explicitee 
a Taide de la figure la sur laquelle les etapes realisees. par le dispositif 
permettant revaluation de la quality des images en sortie d'un processus de 
codage, sont obtenues par differents dispositifs traitant, d'une part les 
signaux d'une image source (10a) et d'autre part les signaux representatife 
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d'une image d6codee (10b). Chaque image est representee par une plurality 
de pixels py arranges en matrice comme repr6sent6s sur la figure 1c et dont 
la taille depend de la definition souhaitee pour I'image. A un pixel py donn6 
correspond une taille de details caract6ristiques exprimee en fonctlon de la 
taille de la matrice definissant le nombre de pixels dans I'image, d'une part 
par la frequence en cycle et d'autre part par la vitesse en degre par seconds. 
Chacune des etapes du proced§ selon la premiere variante de realisation est 
appliqu6e § la fois d I'image source et a I'image decodee. Par Image 
d6cod6e. il faut entendre toute image vid6o obtenue en sortie d'un dispositif 
de codage d6codage permettant la transmission selon un standard tel que, 
par exemple, MPEG. Pour mieux fixer les idees, le lecteur pourra egalement 
se referer ^ I'annexe 1 qui repr6sente les differentes 6tapes num6rot6es de 
1 a 7 du proced6 mis en oeuvre selon la premiere variante. 

"On premieraispositifTlarib) de^pretraitement'des images sources 



15 (10a) et d6cod6es (10b) permet de mettre en ceuvre une premiere etape de 
traitement dite etape de pretraitement r6alisant une correction Gamma des 
signaux repr§sentatife de I'image et une correction du contraste par une loi 
de Weber. La loi de Weber tient compte du fait que roeil est sensible aux 
contrastes et que lorsqu'on regarde un spot lumineux d'intensit6 I + dl sur un 

20 fond ayant une intensite lumineuse I, le rapport dl/l appele rapport de Weber 
est pratiquement constant a environ 2% sur une large plage de seull 
lumineux excepts pour les intensites lumineuses tres faibles et les intensltes 
lumineuses tres elevees. La correction de contraste prend en compte une 
forme de saturation de I'oeil liee au fait que, par exemple, une zone 

25 d'intensite moyenne a cote d'une zone d'intensite forte sera moins bien 
distinguee qu'une zone d'intensite faible A cote d une zone a intensiie 
moyenne. 

Pour tenir compte de cet effet, on corrige le signal de tension qui va 
commander I'intensite lumineuse de I'affichage effectue par le tube 
30 cathodique par une loi dite de Weber exprimee ci-apres : 
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Vout=%^o^iofi+ioo^^ 



_ oiu - T 



dans laquelle Lmctx represente I'intensit6 lumineuse maximum 
approximativement 6gale d cent candela par m§tre carre (X;„ax«100cd/m^) et 
Ldisplay rintensrt6 lumineuse souhalt§e. 

Cette loi mathematique est mise en oeuvre par un dispositif 
electronique permettant d'effectuer ces calculs. A titre d'exemple, un tel 
dispositif peut etre constltu6 d'un microprocesseur assocle a des memoires 
qui contiennent le programme correspondant a I'algorithme de calcul. 

La correction Gamma permet quant a elle de s'affranchir de la 
reponse du tel^viseur, c'est-a-dire des caracteristiques du tube permettant 
j:affichage^.JEn effetJM-dispositifs_d'affich^^ des 
dispositifs non lineaires et I'intenslte de lumi^re reprodulte sur l'6cran d'un 
moniteur cathodique n'est pas, en I'absence de correction, proportionnelle d 
sa tension d'entree. La correction Gamma est un processus de 
compensation de cette non linearite de fagon a obtenir une reproduction 
correcte et proportionnelle de I'intenslte lumineuse correspondante d la 
tension d'entree. L'image d'un ecran est subdivisee en pixels organises en 
matrice dont la position du pixel py est definie par les indices i et j de la 
matrice. La valeur py du pixel est representative de I'intenslte souhaitee. 
Pour corriger les phenom6nes lies au tube cathodique, on applique a cette 
valeur representative de la tension destinee a obtenir I'intenslte lumineuse 
souhaitee, une loi de correction qui con-espond § l'§quation suivante : 



"t 



Gamma (y) ayant une valeur comprise entre 2,3 et 2,6 selon les particularites 
du tube cathodique. Dans cette fonnule. e est le niveau de grille de la valeur 
du pixel Pij, emax est la valeur maximum possible pour e, par exemple 256, si 
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les signaux de commande sont exprim6s sur 8 bits, et Lmax est Tintensite 
correspondant a e^ax en cd/m2, Lmax etant approximativement egal a 
100 cd/ml 

Une autre formulation de la loi de Gamma peut etre la suivante : 
5 y = Ks V 

avecV=kaE^ 

dans laquelle y est la luminance, V la tension de luminance, E I'eclairement 
de rimage analys§e ,y est un exposant voisin de 2,2 pour les tubes images 
en noir et blanc et ya correspond d une valeur de 0,45 admise couramment 

10 pour la t6l§vision couleur, Ks et ka des coefficients de proportionnalite. 

Cette operation de correction Gamma et de Weber permet de 
transformer la valeur du pixel re?ue en entr6e du circuit de pr§traitement (1a) 
en une valeur finale P'ij=lg(Pij) py. Ig qui suit la loi correspondant a la courbe 1 
representee S la figure lb. Chacune des images sources (10a) donne 

15 naissance a une plurality de pixels pr6trait6s p'ag et chacune des images 
d6cod6es (10b) donne egalement naissance a une seconde plurality de 
pixels pr6trait6s p'ty. 

Pairall6lement d cette operation de traitement, une deuxieme etape 
(2a, 2b, fig.1) dite d'estimation de mouvement permet la construction pour 

20 chaque image source et decodee de I'image champ de mouvement, est mise 
en oeuvre a partir de chaque sequence d'image. Cette construction du 
champ de mouvement s'effectue entre t et t-1 par des calculs classiques, tels 
que ceux faisant appel a la methode explicit6e dans le livre publi6 par Don 
Pearson aux editions Mac Graw Hill Book Company et intitul6 "Image 

25 processing" "The Essex series in telecommunication and information 
systems", page 47 et suivantes. Cette estimation de mouvement sur une 
sequence d'images peut utiliser soit la methode diff6rentielle, soit la m6thode 
de coincidence de blocs (Block matching), soit la methode de Fourier ou la 
methode d'estimation en espace tridimensionnel. Pour chaque image, on 

30 obtient ainsi un certain nombre de vecteurs de mouvement (Vj) et I'image 



peut etre ainsi d^coupee en regions (Ri,...,RiVi,...RnVn, fig.ld) basee sur 
I'information mouvement, cheque region (RjVj) est done caracterisee par le 
fait que tous les pixels de cette region ont un seul vecteur vitesse (Vj) 
associe. Cette repartition de I'image en region a vitesse constante constitue 

5 une troisieme etape (3a, 3b) dite de segmentation du champ de mouvement 
appiiqu^e § chaque image source et decodee. Un champ de mouvement est 
done rendu homogene grace a la technique de segmentation et les pixels de 
mouvement proches sont regroupes en une meme region. Ce champ de 
mouvement ainsi segmente est alors plus proche du mouvement vrai des 

10 objets de la sc6ne representee par les images. Cette homogeneisation ^ 
permet ainsi d'effectuer un leger debruitage et d'avoir un nombre r6duit de 
mouvements diff6rents correspondent a un petit nombre de vitesses Vj ce 
qui reduira le nombre de filtres a-calculer et-^-m6moriser-§-r6tape suivante - - — 
pour eviter d'avoir (255)^ mouvements ^ calculer et a memoriser dans le 

15 cadre d'une image constituee de 255*256 pixels Py. L'estimation du 
mouvement peut egalement utiliser la technique fondee sur I'extraction 
d'objets particuliers dans la scene, tels que des bords ou des coins d'objets 
consistent § suivre le mouvement de ces elements particuliers d'uhe image a 
{'autre Ceci foumit des informations de mouvement S diff§rents endroits de 

20 I'image et un processus d'interpolation est utilise pour affecter des vecteurs ^ 
de mouvement aux aires d'images restantes. Une fagon de mesurer les 
mouvements des angles ou bords d'el6ments consiste a appliquer un filtre 
passe haut a I'image pour isoler les bords et ensuite utiliser la technique 
fondee sur la methode des differentiels pour mesurer la valeur du 

25 mouvement. La formation de bordures peut etre attenu6e a I'aide d'un filtre 
pesse bas pour r^duire les effets de bruits et pour permettre la mesure de 
mouvements importants. Un filtre passe bas peut. par exemple, etre 
constitue pour un espace metriciel de dimension 3.3 par la matrice suivante: 
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1 


1 


1 


1 


1 


1 


1 


1 


1 



Un filtre passe haut peut, par exemple, etre constitue pour un 
espace matriciel de dimension 3.3 par la matrlce suivante : 



-0,125 


-0.125 


-0,125 


-0,125 


1 


-0,125 


-0.125 


-0.125 


-0.125 



5 

Ainsi, une table donnant pour une vitesse estim6e v les pixels pg de 
rima ge qui ap partiennent a ce tte region ayant cette Vitesse estim6e v en 
cycle par degre est memorisee sous forme de la matrice de la figure 1d dans 
les moyens de memorisation du dispositif de segmentation du champ de 

10 mouvement en vue de son utilisation ulterieure. 

Le processus se poursuit par une quatrieme etape mise en oeuvre 
en parallele pour chacune des images sources (10a) et decbd^es (10b) 
consistant a construire pour chacune un filtre psychovisuel (4a, 4b). Cette 
quatrieme etape s'effectue, par exemple, a partir d'une base de donnees de 

15 courbes de filtre comportant une pluralite de courbes telles que celles 
representees a la figure 3 exprimant Tinfluence du facteur humain H en 
fonction de la vitesse du mouvement Vj et de la frequence spatiale fj, A ces 
valeurs Vj. fj correspond une valeur de filtrage H. Si la vitesse Vj est comprise 
entre les vitesses v^ et Vg de deux courbes (Hvi.Hv^), le dispositif effectue 

20 une interpolation pour determiner la valeur H correspondante. Cette 
interpolation peut etre, par exemple. lineaire. Cette etape peut egalement 
s'effectuer par calcul direct de la valeur H pour une vitesse donn6e ^ partir 
d'un modele analytique de Tinfluence psychovisuelle. Le modele analytique 
du filtre peut. par exemple, etre represents par la formule suivante : 




G(a,v) = [6.1 + 7.3| log(v/3)l ^] x va^ exp[-2a(v+2)/45.9] avec a=27rf. 

A la cinqui^me etape, le dispositif synthetise une valeur de filtrage 
h(s,v) dans le domaine spatial pour chaque valeur de filtrage H(fs.v) dans le 
domaine frequentiel associee ^ une vitesse v et une frequence % 
appliquant une transfornnee rapide de Fourier inverse (FFT^) sur les valeurs 
de filtrage H(fs,v) dans le domaine frequentiel, ce qui s'exprime par 
i'expression h(s,v) = FFT^IH(fs.v)] oil fs repr§sente la taille caract^ristique 
des details en cycle par degr6, v leur mouvement exprime en duree par 
seconde et s la resolution. 

La valeur h du domaine spatial ainsi d6terminee pour chaque image 
source (5a) et d6codee (5b) dans la cinquieme etape sera appliqu6e au 
cours d'une sixieme etape (6a, respectivement 6b) d chacun des pixels 
pr6traites Pay de I'image source sortant du pr^traitement et de I'image 
decodee p'bg sortant du pr^traitement. 

De cette sixieme etape (6a, 6b) de filtrage, il r6sulte d chaque fois un 
couple de valeurs de pixels filtr6es p'a,y et p'tfy qui est ensuite utilise dans 
une derniere etape (7) de construction d'une matrice des disparites (Djsp) en 
effectuant le calcul de I'erreur quadratique entre chaque couple de valeurs 
de pixels. 

Dtep = (p'afij-P'bfii)" avec n = 2 ou autre 
Cette matrice donnera ainsI une appreciation objective des 
distorsions perceptibles a I'oeil humain introduites par le processus de 
codage decodage que Ton veut estimer. 

La deuxieme variante de realisation de I'invention va maintenant etre 
explicit6e en liaison avec la figure 2 qui represente les differents moyens 
permettant la mise en oeuvre du proc6d§ selon cette deuxieme variante. 
Pour mieux fixer les id6es, le lecteur pourra egalement se referer a I'annexe 
2 qui represente les differentes etapes numerot§es de 1 a 1 1 du proced6 
mis en oeuvre selon la seconde variante et qui represente egalement la 
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diminution de taille obtenue par les operations de construction des 
pyramides de Laplace. 

Dans cette variante, les quatre premieres etapes de la premiere 
variante sont appliquees a Timage source (10a). a savoir la pr6correction 

5 (1a), la construction de champs de mouvement (2a) et la segmentation (3a) 
ainsi que la construction du filtre psychovisuel (4a), Les pixels Pag provenant 
du traitement de I'image source (10a) et resultant de Tetape de pretraitement 
et les pixels Pby resultant de Tetape de pretraitement de Timage decodee 
(10b) sont chacun soumis dans une etape (5.1a et respectivement 5.1b) a 

10 un filtrage de decimation Vz (F1/2). Ce filtrage est un filtrage passe bas qui 
permet a partir d'une matrice de pixels representant une image Pn.id'un 
niveau n-1 donn6 d'obtenir Timage Pn de niveau n suivant. Ceci s'exprime 
par la relation : 

Pn = F1/2 (Pn.i) avec n > 0 

15 Po etant I'image originale. 

A titre d'exemple le filtre de decimation peut etre constitu6 pour un 
espace matriciel 3*3 par la matrice suivante : 



1 


2 


1 


2 


4 


2 


1 


2 


1 



20 Cette operation de filtrage par decimation % par un dispositif de 

calcul a pour r§sultat de r6duire une matrice de pixels de taille m-n 
representant I'image source Pgo en une matrice Pgi de niveau 1 de taille 
m/2.n/2, la matrice Ps„ de niveau n etant de taille m/2n-n/2n. De menne, cette 
operation de filtrage par decimation !4 a pour resultat de reduire une matrice 

25 de pixels de taille m.n representant I'image d6cod6e ou con-ompue Pjo en 
une matrice P^i de niveau 1 de taille m/2-n/2, la matrice Pjm de niveau n 
etant de taille m/2n.n/2n. Le dispositif de calcul memorise done pour chaque 



image source et d6codee le niveau n et le niveau n+1 suivant dans sa 
m^moire. 

Ensuite, d l'§tape suivante {5.2a, respectivement 5.2b), le dispositif 
de calcul retranche a chaque image de niveau n, I'image Pn+i, de niveau 
immediatement suivant, expans^e par 2 (E2) pour obtenir ce qu'il est 
convenu d'appeler une succession de matrices constituant ce qu'il est 
convenu d'appeler une pyramide de Laplace Ln, selon la formule : 

Ln = Pn - E2 (P„+i) pour n < N 

avec Ln= Pn 

L'operation d'expansion effectu6e par E2 consiste S interpoler I'image 
Pn+i (de taille m/2*n/2) pour obtenir une image de taille m*n. 

Cette operation d'expansion ou d'interpolation met en jeu plusieurs 
-matriees-d'interpolation-en fonction de la position du pixels interpoler-. - - 

Ceci permet de constituer, pour I'image source, une pyramide LS„ de 
matrices m§moris6es et pour I'image decodee ou corrompue une seconde 
pyramide LD„ de matrices m6moris6es. Selon le choix du filtre fA , I'image 
LSn obtenue d la fin de l'6tape ci-dessus est une bonne approximation de 
r^nergie comprise dans une bande de frequence centree autour de la valeur 
fp = 1/(n+1). Pour plus de details sur les pyramides de Laplace ou encore 
dites gaussiennes et sur les matrices d'expansion, le lecteur peut se reporter 
a I'article " the Laplacian pyramid as a compact code image " paru dans la 
revue IEEE transactions on communications VOL. COM. 31, N" 4, Avril 1983 
pages 532 a 540 auteurs P. J. Burt et Ed. H. Adelson. 

Les pyramides de Laplace (LSp, LDJ sont obtenues par les etapes 

(5.2a. et 5.2.b). 

A I'etape suivante, le dispositif constrult, S partir des regions (RjVj) 
caracterisees par le m§me vecteur vitesse (Vj) et definies en sortie de i'etape 
de segmentation (3) du cliamp de mouvement appliquee d I'image source, 
une pyramide multir§solutions Rn de I'image region en partant de I'image 
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region originale Rq et en appliquant a cette image region originale Rq un filtre 
median Gvj de decimation Vz selon la formule : 



Ro = i'image region originale. 

Pour plus de renseignement sur la constitution d'un filtre median 
connu de I'homme de m6tier, le lecteur pourra se referer au chapitre 4 du 
livre intitule " non linear digital filters, principles and applications " publie en 
1990 aux Editions Kluwer Academic Publishers par I. Pitas et A. N. 
Venetsapoulos. 

La pyramide de mouvement (multiresolution) est r^aiis^e 
uniquement pour I'image source (10a) et la valeur du pixel de LS„ et LDn 
represente I'energie qu'il y a dans une bande de frequence alors que par Rp 
on a le mouvement. 

bans cette etape, le dispositif de calcul n'effectue pas le calcul 
correspondant d I'etape 5.2, c'est-a-dire a la soustraction de I'image Rn 
correspondant a chaque niveau n de I'image Rn+i du niveau n+1 
imm6diatement suivant expans§ par 2. A cette 6tape de calcul, les matrices 
constituant la pyramide multiresolution sont m6moris6es et permettent 
d'obtenir pour chaque pixel de chaque niveau n de Rp une valeur locale du 
mouvement. Le filtre m6dian d'une fen§tre de n.m pixels est obtenu en 
classant les valeurs des pixels et en retenant le pixel ayant la valeur 
m^diane. 

L'application du filtre median Gvj ^ la matrice de pixels Rp de 
dimension m*n permet d'obtenir une matrice Rn+i de dimension m/2*n/2. 
L'operation de decimation est incluse dans celle du filtre median. Le filtre 
median a le meme effet sur les images Rn que le filtre F% sur les images P„ : 
il r6duit leur taille par 2 en horizontal et en vertical, sauf qu'au lieu d'etre un 
filtre matriciel classique, c'est un filtre "median", c'est-a-dire bas6 sur une 
analyse de la statistique locale. 



Rn ~ G (Rn-l) 



Les pyramides de Laplace (LSn, LDn) calcul6es a I'etape 5.2 sont 
ensuite utilis6es dans une etape de calcul des differences inter pyramides 
niveau ^ niveau representee a I'etape 7 salon la formule : 

Diff„=LSn-LDn. 

5 Ceci permet d'obtenir des matrices Diffp dont chaque coefficient 

exprime les valeurs des differences des coefficients des matrices des 
pyramides de Laplace source (LSJ et decodee (LDp) pour le meme niveau n 
et ceci pour chaque niveau de 0 ^ n. Dans la pyramide de Laplace LSn 'a 
valeur du pixel repr6sente I'^nergie qu'il y a dans une bande de frequence. 

10 En faisant la difference inter pyramides LSn - LDn on a pour une bande de 
frequence donn§e la disparity fr6quentielle entre les deux images. 

Ce r6sultat est en fait pondere par la sensibility de I'oeil pour cette 
-frequence-qui-est-exprim6e par la-felation-de-l-influenee relative de ractlvite 
de la frequence fn. Cette influence relative de I'activite de la frequence fn peut 

15 etre masquee par une activite importante dans les frequences superieures. 
Pour determiner et tenir compte de cette influence relative de I'activite dans 
une etape de masquage (8), le dispositif de calcul commence par 6valuer 
I'influence locale £„ d'un pixel pg qui est definie par la valfeur du resultat du 
calcul des differences inter pyramides applique au pixel pg, ce resultat etant 

20 eieve §1 la puissance q : ^ 

^«=(Diffn(Pij))^ 

cette valeur d'influence locale permet au dispositif de calcul de determiner 
une matrice exprimant I'influence relative de I'activite par un circuit qui met 
en oeuvre la formule suivante : 

^ " Y.^{E,y ~ " 

k<n 

avec m(Ek)= Ek si Ek > S 
et m(Ek) = S si Ek < S 

avec par exemple S=0,5% (valeur maximale possible de E^) 




Comme dans la variante precedente, le dispositif de calcul effectue 
une etape de filtrage 4 par Tutilisation d*une base de donnees (BD) 
contenant une pluralite de courbes de filtrage exprimant {'influence du facteur 
humain sur la perception visuelle des images. Ces courbes de filtrage 

5 permettent, a partir des valeurs de la frequence et de vitesse correspondant 
a un pixel py de determiner un coefficient H de ponderation pour ce pixel. 
Ainsi, pour chaque pixel py de la matrice correspondant a une region Rn 
de Vitesse, le dispositif de calcul determine une valeur H qui vient ponderer 
rinfluence relative In. Cette etape de ponderation (9) est obtenue par un 

10 dispositif de calcul mettant en ceuvre Tequation : 

La mise en ceuvre de cette equation permet d'obtenir une pyramide 
de matrices. Quand le programme du dispositif de calcul selectionne un pixel 
dans une matrice de niveau n de la pyramide de Laplace, d ce niveau n 

15 correspond une frequence spatiale fp et le dispositif de calcul salt associer le 
pixel de Timage Rn auquel correspond une valeur de vitesse v. Par utilisation 
des bases de donnees et des courbes enregistrees dans ces bases de 
donnees, le dispositif de calcul determine soit directement, soit par 
interpolation entre deux courbes, la valeur du coefficient de gain H. dans 

20 cette deuxleme methode, on travaille directement sur des objets (pyramides 
de Laplace) qui correspondent a des grandeurs frequentielles. On n'a done 
pas besoin de passer dans le domaine spatial, car tous les calculs se font 
dans le domaine frequentiel. Cette valeur H va venir ponderer rinfluence 
relative de Tactivite (In). Cette etape de ponderation (9) permet d'obtenir une 

25 pyramide de matrices auxquelles peut s'appliquer une etape (10) facultative 
assurant le filtrage directionnel pour tenir compte des directions 
psychovisuelles privll§giees par le regard humain. 

Ainsi, on peut filtrer les Images constitutes par les matrices Tn par 
des filtres directionnels favorlsant une direction par rapport ^ d'autres, ces 

30 filtres etant constitues de matrices de coefficients de dimension n.n 



correspondant d la dimension du facteur d'images T„. Un exemple de 
matrice de filtre directionnel a 0' est donne ci-apres pour une dimension 5/5. 



0 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0.094 


0,141 


0,28 


0,141 


0,094 


0 


0 


0 


0 


0 


0 


0 


0 


0 


0 



Un exemple de matrice de filtre directionnel a 90"* est donne ci- 



apres. 



0 


0 


0,125 


0 


0 


0 


0 


0.188 


0 


0 


0 


0 


0,375 


0 


0 


0 


0 


0,188 


0 


0 


0 


0 


0,125 


0 


0 



Un exemple de matrice de filtre directionnel a 45*" est donne ci- 



10 



apres. 





0,016 


0 


0 


0 


0 




0 


0,023 


0 


0 


0 


0 


0 


0,047 


0 


0 


0 


0 


0 


0,023 


0 


a 


a 


a 


0- 


0t046- 

















Le resultat de cette 6tape (10) de filtrage directionnel est envoye sur 
un circuit sommateur pour recomposer dans une etape (11) les pyramides 
15 multiresolutions P'n par l'§quation 
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PV = E2 (PVi) + T„ (n<N) 

avec P N = Tn 

Ceci entraTne done 

Pn.i = E2 (Tn)+Tn.i 

5 Le processus est reiter6 jusqu'a obtenir P 0 qui represente la matrice 

constituant la carte des disparites. 

Ces etapes des deux variantes de realisation de Tinvention sont 
mises en oeuvre a Taide de circuits d microprocesseurs ex6cutant les 
programmes appropri6s permettant de r6aliser les etapes exposees 

10 pr6cedemment. Ces circuits comportent egalement des moyens de 
memorisation pour memoriser en plus des programmes a executer, les 
matrices de pixels ou les matrices de regions ou encore les resultats 
intermediaires permettant d'attendre Tetape de calcul suivante ou encore les 
filtres ^^ppllquer"aux re intenn6diaires ouTfnaux! 

15 Ceci est utilise pour comparer les performances subjectives de 

differents appareils codeurs ou encore pour comparer la performance 
subjective de differents algorithmes de codage et/ou effectuer la mesure de 
la perception d'artefacts dus a des traitements d'image. En fonction de ces 
performances, le circuit de calcul mettant en oeuvre une des deux variantes 

20 de realisation de Tinvention peut modifier le codage en faisant. par exemple. 
retroagir la mesure de quaiite subjective des images codees ainsi effectuees 
sur les parametres globaux et/ou locaux du codage. Un parametre global sur 
lequel s'effectue la retroaction peut. par exemple, etre le debit moyen et un 
parametre local peut, par exemple, etre le pas de quantification local utilise 

25 lors du codage. Cette retroaction peut etre dynamique en cours de codage. 
les erreurs retroagissant sur le pas de quantification local, la taille des 
images, la forme du GOP (Group of pictures) etc. La retroaction peut 
egalement s'effectuer de fagon iterative dans le cas des codages pour 
disque video (DVD) ou CD-ROM. Dans ce cas, tant que I'erreur est 

30 superieure ^ un seuil et/ou n'est pas homog^ne sur toute Timage, la 




retroaction du circuit de calcul commande une diminution de la severite des 
parametres de codage globalement et/ou localement. Tant que Terreur est 
inferieure a un seuil, la mesure de la qualite subjective des images 
effectuees par le circuit de calcul pemnet d'augmenter la severite des 
parametres de codage. 

Enfin. une avant-derniere utilisation de la methode d'estimation ou 
de mesure de la qualite subjective d'une sequence d'images issues d'un 
codeur peut concerner le codage oriente objet. Dans ce cas, Tutilisation du 
dispositif de calcul de Tune des methodes mise en oeuvre par ce dispositif 
permet d'assurer une qualite subjective constante des divers objets de la 
scene ou d'assurer une qualite subjective relative inter-objets donnee. Enfin, 
le precede et le dispositif d'estimation peuvent permettre de modifier les 
conditions de reallocation dynamique des debits alloues a chacun des 
canaux d'une diffusion avec multiplexage statistique pour assurer une qualite 
subjective donn6e et homogene des programmes diffuses. 

D'autres modifications a la portee de Thomme de metier font 
egalement partie de Tesprit de I'invention. 
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ANNEXE 1 



Methode 1 

5 

Etape 1 Pre-Correction des images: Correction Gamma 
de Tecran et Correction du Contrast (Loi de Weber). 

Etape 2 Construction de 1' image Champ de Mouvement 
10 a partir de la Sequence Source et pour chaque image. 

Etape 3 Segmentation du Champ de Mouvement. Pour 
chaque image, on dispose ainsi d'xme Segmentation en Regions basee 
sur r information Mouvement. 
15 Chaque Region(v) est done caracterisee par un vectexir 

Vitesse v. 

Chaque pixel de chaque image (Source ou Decodee) 
appartient a une Region correspondant a une vitesse estimee v (en 
cycles par degre). 

20 

Etape 4 Pour chaque Region(v), Construction du Filtre 
Psychovisuel correspondant, a partir d'lme BDD Filtres { H(fs 
,V2)i=^]^„.N } interpolation des filtres. 

25 Etape 5 Pour chaque Region(v), Synthese du Filtre en 

Spatial par FFT inverse: 

h(s,v) ^ FFTi [H(fs ,v)] 

Etape 6 Filtrage des images Soxirces et Decodees pour 
30 obtenir deux autres images: 




SourceF et DecodeeF. 

Chaque pixel P de I'image Sotirce/Decodee est filtre par 
le Filtre h(s,v), correspondant a la Region(v) a laquelle appartient P, 
centre sur P et applique sur Timage Source/Decodee. 



Etape? Construction de la Carte des Disparites ou 
Erreurs psychovisuelles. 

Err = (SourceF- DecodeeF)^ (n^2, ou autre) 
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ANNEXE 2 



M ahode 2 



10 



Etape 1 (voir Methode 7). 

Etape 2 (voir Methode 1). 

Etape 3 (voir Methode 1), 

Etape 4 (Vo/r Methode i). 



15 Etape 5 Decomposition des images Sources et Decodee 

en Pyr amide de Laplace de N niveaux construite en dexxx etapes: 

• Chaque niveaux est d'abord obtenu par Filtrage- 
Decimation Vt, (Passe-bas) du niveau immediatement precedent 
(Pyramide multi-resolution, figure 4). 

20 

Pn = F^/P„.0 n>0 
Pq = image originale 

• Puis on retranche a chaque niveau Pn le niveaux 
25 immediatement suivant expanse par 2 pour obtenir Ln (Pyramide 

de Laplace, figure 5). 



L„ = Pn-E2(P„+i) n<N 
L^^=PN 
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Ce calcxil permet d'obtenir une representation d'xine pyramide multi- 
resolution PnConforme a la figure 4 et-une representation d'une pyramide de Laplace 
LnComforme a la figure 5. 

5 Si le Filtre Fy, est bien choisi, Timage L,, est une bonne 

approximation de Tenergie comprise dans une bande de frequences 
centree autour de 
f =l/(n+l). 

Finalement, on dispose de deux Pyramides de Laplace: LS„ 
10 (Source) et LDn (Decodee). 

Etape6 Sur le meme principe, on construit la Pyramide multi- 
- -resoliition-- Rn--de-J:image_RegionSo 
Decimation V%, Ainsi, pour chaque pixel de chaque niveau des Pyramides de Laplace 
15 (Etape 5), on dispose de la valeur locale du mouvement. 

Etape 7 Calcul des Differences inter-Pyramides niveau a niveau: 

Diffn-=LSn'LDn 

20 Etape 8 Application du principe de masquage frequentiel (Texture- 

Masking^: 

L'influence relative de Tactivite a la frequence f„ est masquee 
par une activite importante dans les frequences superieures (fk<n)- 

On definit alors Tinfluence locale relative du pixel pi, In(pi) 

25 par : 

E„ avec: Ey^ 

avec q==2 par exemple. 



Etape 9 Filtrage des Pyramides de Laplace Source et Decodees: 
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Chaque pixel pi de Ln est pondere par la valevur H(fn,v), 
correspondant a la Region(v) a laquelle appartient pi dans R^, et par 
rinfluence relative In. 

T„(pi) = I„(pi)xH 

Etape 10 Filtrage Directionnel: pour tenir compte des directions 
psychovisuelles privilegiees par ie regard humain, on peut filtrer les images Tn par 
des filtres directionnels favorisant une direction par rapport a d'autres. 

Etape 1 1 Construction de la Carte des Disparites ou 
Erreurs psychovisuelles: On recompose les pyramides multi- 
resolution 

P'n = E2(PW + T„ (n<N) 



La Carte des disparites conrespond a P'o • 



REVENDICATIONS 

1. Proc§d6 d'evaluation de la quality d'images codees caracterise en 
ce qu'il comporte : 

5 a) une etape de traitement du signal repr6sentatlf de I'lmage pour 

obtenir un signal traite, 

b) une §tape de construction a partir du signal representatif de 

I'image cod6e, d'un signal representatif de I'image champ de mouvement a 

partir de la sequence source, 
10 c) une §tape de constitution d'un signal representatif de la 0 

segmentation du champ de mouvement et de stockage des pixels d'images 

repr6sentatifs de chaque region ayant un champ de mouvement different a 
-une-adresse definie-par-rapport aux vecteurs vitesse-estim6s-dans I'etape de 

construction du champ de mouvement permettant de d6temiiner les pixels 
15 ayant des vecteurs Vitesse differents, 

d) une 6tape de d6temnlnation ou de calcul d'un filtre humain 
psychovisuel ^ appliquer en fonction de la Vitesse estim§e de la region. 

e) une etape de filtrage du signal traite, et 

f) une etape de construction de la carte des disparit6s entre les 

20 signaux representatifs de I'image obtenus apres i'etape de filtrage et les ^ 
signaux representatifs de I'image d§cod6e obtenus apres I'etape de filtrage. 

2. Precede d'evaluation de la qualite d'images codees selon la 
revendication 1. caracterise en ce qu'il comporte une etape conslstant a 
appliquer chacune des etapes precedentes a I'image source et ^ I'image 

25 decodee. 

3. Precede d'evaluation de la quality d'images codees selon la 
revendication 1, caracterise en ce qu'il comporte une etape de 

decomposition frequentlelle des images (FFT, sous-bande ) precedant 

retape de filtrage constitu6e d'une ponderation par un coefficient deduit de 

30 courbes prenant en compte la Vitesse estimee et la bande de frequence 




consideree, pour tenir compte de I'influence relative de la vitesse et de la 
frequence spatiale sur la perception des images animees. 

4. Proc6d§ selon la revendication 1 , caracterise en ce que I'etape de 
filtrage psychovisuel s'applique a des matrices representatives des 

5 differences inter-pyramides entre les pyramides de Laplace des images 
sources traitees et celles des images d6cod6es trait§es apr§s ponderation 
par. d'une part I'influence locale representative de la frequence du pixel 
concerne et d'autre part, un coefficient de filtrage deduit de courbes de 
filtrage prenant en compte la vitesse estimee et la bande de frequence 

10 correspondant au niveau de la pyramide de Laplace a laquelle le pixel 
appartlent dans une pyramide multir6solution obtenue par construction d'une 
pyramide ^ partir de I'image de chaque region de vitesse diff6rente. 

5. Procede selon la revendication 1, caracterise en ce que les 
cdurbeVde filtrage psychovisl^^ par une succession de 

15 courbes arrang6es sous fomrie de base de donnees et memorisees dans le 
systeme. et eventuellement interpolation ^ partir de ces courbes, soit 
obtenues par une representation analytique mise en oeuvre par des moyens 
de calcul permettant de calculer chaque courbe. 

6. Procede selon la revendication 4, caract6ris§ en ce que I'fetape de 
20 construction de la carte des disparites s'effectue par recomposition des 

pyramides multiresolutions filtr^es obtenues a I'etape precedente. 

7. Procede selon la revendication 4 ou 6, caracteris6 en ce que 
I'etape de traitement de I'image comporte une etape de decomposition des 
images sources et decodees en pyramide de Laplace de n niveaux et une 

25 etape de construction de la difference Inter-pyramide. 

8. Procede selon la revendication 1 , caracterise en ce que la vitesse 
ou valeur locale du mouvement est obtenue par une construction eventuelle 
de filtres puis application du filtre construit ou par application d'un filtre 
median. 




9. Procide selon la revendication 1 ou 4, caracterise en ce qu'il 
comporte une 6tape de precorrection des images en effectuant une 
correction Gamma et une correction par la loi de Weber. 

10. Precede selon la revendication 7, caracterise en ce que la 
5 correction par Gamma y est la suivante : 

y= KsV'^avecV=kaE^^ 
dans laquelle y est la luminance, V ia tension de luminance. E Tfeclairement 
de I'image analysee eclairement, ys est une exposant voisin de 2,2 pour les 
tubes images en noir et blanc et a une valeur de 0,45 admise couramment 
10 pour la television couleur. 

11. Proc6d6 selon la revendication 1, caracteris6 en ce que le 
filtrage est obtenu par construction du filtre psychovisuel correspondant a ia 
Vitesse estimee a partir d'une base de donnees de filtres et interpolation 
entre les deux filtres con-espondant aux regions les plus proches de la region 

15 dont la Vitesse a 6t6 estim6e. 

12. Precede selon la revendication 4, caracterise en ce que 
I'infiuence (l„) locale relative du pixel Pi concerne est obtenue par calcul 
d'une valeur Bp representant la puissance q'*""^ de la difference niveau a 
niveau inter-pyramides entre les pyramides sources et pyramides decodees 

20 de m§me niveau du pixel concern6. 

13. Precede selon la revendication 12, caracterise en ce que le 
calcul de !„ est effectue par utilisation de la fonnule suivante : 

k<n 

av e cC„^(Diffn(p,j))\ 

25 m(Ek)= Ek si Ek > S 
et m(Ek) = S si Ek < S 

avec par exemple S=0,5% (valeur maximale possible de Ek) 
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14. Procede selon la revendication 4, caracterise en ce que le 
filtrage comporte un filtrage directionnel des images selon une direction 
determinee plutot qu'une autre. 

15. Procede selon la revendication 9, caract6ris6 en ce que la 
correction Gamma s'effectue par un dispositif de calcul mettant en oeuvre 
{'Equation suivante : 



^display ^max 



e 



e etant la valeur de niveau de grille du pixel, e^ax ^tant la valeur maximum 
exemple 256 si le codage s'effectue sur 8 bits, L^ax 6tant Tintensit^ 
10 correspondant a e^ax cd/m . 

16. Proced6 selon la revendication 9. caracteris6 en ce que la loi de 
Weber est mise en oeuvre par un dispositif de calcul qui realise la fonction 
suivante : 



15 17. Precede selon la revendication 1 ou 4, caract6rise en ce que le 

calcul du filtre s'obtient par la formule suivante : 

G(a,v) = [6.1 + 7.3| log(v/3)| ^] x va^exp[-2a(v+2)/45.9] 
avec a=27if, f=frequence spatiale, v=vitesse. 

18. Utilisation du proc6de selon une des revendications precedentes 
20 dans un dispositif de codage caracterisee par une retroaction dynamique en 

fonction des disparites psychovisuelles calculees par le dispositif de calcul 
mettant en ceuvre le proced6 sur Tun des parametres utilises par le dispositif 
de codage au cours du codage. 

19. Utilisation du procede selon la revendication 18. caracterisee en 
25 ce que les disparites calculees sont comparees a un seuil pour modifier les 

parametres de codage de Tappareil de codage jusqu'^ ce que le seuil 
souhaite soit franchi. 



20. Utilisation du proced6 selon la revendication 19, caracteris6e en 
ce que I'un des parametres est soit le pas de quantification, solt la taille des 
images, soit la forme du groupe d'images GOP. 

21. Utilisation du proc6d6 selon la revendication 18, caracteris6e en 
ce que i'homog6neit§ des disparites calculees est analysee par le dispositif 
de calcul pour agir sur les parametres de codage de I'appareil de codage. 

22. Utilisation du proced6 selon les revendications pr6c6dentes, 
caracterisee en ce que les parametres de codage des objets differents d'une 
Image dont le codage est oriente objet sont modifies en fonction d'une 
disparite souhait6e constante. 

23. Utilisation du proced6 selon les revendications 18 22, 
caracterls6e en ce qu'elle consiste S effectuer une reallocation dynamique 
des d6bits-alloues-a=un appareil-de-GGdage-aveG-muttiplexager— 

24. Dispositif d'6vaIuation de la qualite d'images cod6es caract§ris§ 

en ce qu'il comporte : 

- un moyen (la, lb) de traitement du signal representatif de I'image 
source (10a) et de I'image d6cod6e (10b) pour obtenir un signal d'image 
source trait§e et d'image decodee traitee, 

- des moyens (2a, 2b) de construction a partir du signal representatif 
de chacune des images, d'un signal representatif de I'estimation du champ 
de mouvement a partir de chacune des images des sequences sources et 
d^cod^es, 

- des moyens (3a, 3b) de constitution d'un signal representatif de la 
segmentation du champ de mouvement et de stockage des pixels d'images 
representatifs de chaque region Rj ayant un champ de mouvement different 
d une adresse d§finie par rapport aux vecteurs vitesse estim6e dans I'etape 
de construction du champ de mouvement permettant de determiner pour 
chacune des images sources et d6codees celles ayant des vecteurs vitesse 
differents, 
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- un moyen (4,5) de determination ou de calcul d'un filtre hunnain 
psychovisuel a appliquer en fonction de la vitesse estimee de la region, 

- des moyens (6a. 6b) de filtrage appliques a chacune des images 
sources traitees et d6codees traitees et 

5 - un moyen (7) de construction de la carte des disparites entre les 

signaux repr6sentatifs de I'image source traitee obtenus apres r6tape de 
filtrage et les signaux representatifs de Timage decodee traitee obtenus 
apr^s l'6tape de filtrage. 

25. Dispositif selon la revendication 24, caracterise en ce que les 
10 moyens de filtrage psychovisuel s'appliquent a des matrices representatives 

des differences inter-pyramides calculees par des moyens de calcul entre les 
pyramldes de Laplace des images sources traitees et celles des images 
decodees traitees apres ponderation par d'une part influence locale 
~7epTesentatlS^lie~la^ tin 

15 coefficient de filtrage deduit de courbes de filtrage memorisees ou calculees 
et prenant en compte la vitesse estimee et la bande de frequence 
correspondant au niveau de la pyramide de Laplace a laquelle le pixel 
appartient dans une pyramide multi resolution obtenue par des moyens de 
construction de cette pyramide multiresolution a partir de I'image de chaque 

20 region de vitesse differente. 

26. Dispositif selon la revendication 24. caracterise en ce que les 
moyens de construction de la carte des disparites effectuent une 
recomposition des pyramides multiresolutions filtrees. 

27. Dispositif selon une des revendications 24 a 26. caracterise en 
25 ce que les moyens de traitement. les moyens de constitution, les moyens de 

determination, les moyens de construction, les moyens de filtrage sont 
constitues d'au moihs un microprocesseur associe a des memoires 
suffisantes pour contenir les programmes permettant de realiser les 
differents moyens et pour contenir les bases de donnees et les informations 




intermediaires necessaires au calcul et d I'obtention de la carte 
disparites. 
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